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Ph, Abstract 

The large observed mixing angle in atmospheric neutrinos, cou- 
■ pled with Grand Unification, motivates the search for a large mixing 

between right-handed strange and bottom squarks. Such mixing does 
not appear in the standard CKM phenomenology, but may induce 
significant b ^ s transitions through gluino diagrams. Working in 
the mass eigenbasis, we show quantitatively that an effect on 

CP violation in 4>Ks is possible due to a large mixing between 

sr and bpi, while still satisfying constraints from b — 57. We also 
include the effect of bL-bn mixing proportional to mhfitan fi. In the 
case where mfe/xtan/3 <^ ^SUSY there may be a large effect in Bs 
mixing correlated with a large effect in B'^ — > (pKs, typically yielding 
an unambiguous signal of new physics at Tevatron Run II. 
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1 Introduction 



Flavor physics has seen tremendous progress in the past few years. The 
discovery of neutrino oscillations by the SuperKamiokande [1 , SNO |2] , and 
KamLAND [3j experiments clearly marks a historic event, while CP violation 
has recently been found in two new manifestations: direct CP violation in 
the neutral kaon system [4^ and indirect CP violation in the the B'^ system 
jS]. On the other hand, we still lack insight into the origin of flavor and the 
patterns of masses and mixings. We need to look for any possible hints of 
physics that give us additional insight into these questions. 

One of the major surprises in neutrino physics was the observation of 
(two) large angles. Unlike in the quark sector where all mixing angles in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix are small, both atmospheric 
and solar neutrino oscillations require large angles. An important question is 
whether the presence of large angles will give us new insight into the origin 
of flavor, masses, and mixings. 

It was pointed out in Ref. jB] that the large angles in the neutrino sector 
may imply large angles in the mixing among right-handed down-type quarks 
if they are grand-unified with lepton doublets. Indeed, some 5*0(10) models 
with Pati-Salam type unification of Yukawa matrices suggest that the large 
mixing angles in neutrinos arise from the charged lepton mass matrices, and 
thus also appear in the down-quark mass matrices. In these models, one 
assumes that these new large mixing angles do not appear in the CKM ma- 
trix because the right-handed charged- current interaction is broken at the 
Pati-Salam unification scale. However, the imprint of the large atmospheric 
neutrino mixing angle may appear in the squark mass matrices as a large 
bR-SR mixing effect though radiative corrections due to the large top Yukawa 
coupling. The large solar neutrino mixing angle, however, does not cause a 
significant effect because of the smaller Yukawa coupling for lower genera- 
tions. The new bn-SR mixing in turn feeds into new effects in 5-physics. In 
particular, there may be large new CP-violating effects in b —>■ s transitions 
and enhanced Bg mixing. It has already been noted that CP violation in 
B^ —>■ (f)Ks is a good place to look for new physics effects [7| |HI. 

The time-dependent asymmetry in B^{B^^) 4>Ks was reported recently 
by both BaBar and BELLE. Their measurements differ from the value in 
the J/ipKs final state by 0{1). The standard model predicts that these 
two channels should give the same value. The significance of the difference is 
2.7cr if the measurements from both collaborations are combined- the current 
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world average for sin 2/3 in the J/ipKs channel is 0.734 ± 0.054, while 

in the (pKs channel S^px = —0.39 ± 0.41 0. This report has already 

sparked many speculations JU] • It is not clear if this is a temporary anomaly 
or a genuine new effect. Nonetheless it is important to study how large the 
new CP violation in B^ —>■ (pKs can be and how it is correlated to Bs mixing 
which will be studied soon at Tevatron Run II. 

In this paper, we investigate the size of CP violation in B^ <pKs as well 
as Bs mixing from a potentially large bn-SR mixing. There have also been 
several investigations of B^ — > (pKs within the context of supersymmetry 
(SUSY) [H1III1II2]. Of the above, only Ref. [12] investigated the correlation 
between the measurement of sin 2/3 in B^ — >■ (pKs and Bg mixing. However 
it uses the mass insertion formalism, which is not necessarily appropriate for 
the large mixing that we will consider. In addition, it appeared before the 
recent experimental results, and so it did not seek to reproduce such a large 
shift in sin 2(3. We perform a calculation in the mass eigenbasis, with a goal 
of determining whether supersymmetry can accommodate the central value 
of the recent experimental results for sin 2/3 in B^ — > (pKs- We then explore 
the consequences for Bg mixing. 

We also emphasize contributions to sin 2(3 that arise from a combination 
of (mfe/i tan /3) and bR-SR mixing. These contributions, which we find to 
be important over a wide region of parameter space, are not easily analyzed 
in the mass insertion approximation. Analogous combinations were studied 
in the kaon system ^3]; but to our knowledge these contributions have not 
been thoroughly analyzed with regard to new physics in the i?^ system. 

The outline of the paper is as follows. In the next section we introduce the 
effective field theory formalism for b decay and work out b ^ s transitions. In 
Section 3 we discuss Bg mixing from large 6_r-S_r mixing. Section 4 is devoted 
to the discussion of correlations between the b —>■ s transition and Bg mixing. 
We conclude in Section 5. Details of some calculations are presented in the 
appendices. In Appendix A we show the loop functions, while the hadronic 
matrix elements are estimated in Appendix B. 

2 CP Violation in 6 ^ s Transition 

In this section we briefly review the well-known effective field theory formal- 
ism for 5-physics (for a comprehensive review see which we use to 
calculate the contribution of supersymmetric particles to 6 — > s transitions. 
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Using this machinery we discuss the contribution of a large mixing between 
right-handed squarks to the CP- violating parameter 5'<^x. We also use this 
formalism to address the constraints on the SUSY contribution that come 
from the 6 — > 57 radiative decay. 



2.1 Effective Hamiltonian 

The b ^ s transitions of interest can be described by the following AB = 1 
effective Hamiltonian: 

^eff = E {Qo, + c[o[) + c,o, + + c,o, + (i) 



1=1 



where 



01 = isa''PLCj)iCj-f^PLh) (2) 

02 = isYPLc)ic^>.PLb) (3) 

03 = isYPLb)is^^PLs) (4) 

04 = {sa^PLbj){s,^^PLS^) (5) 
a = {s^^PMs^^Prs) (6) 

a = {s^rPLh){s,-i^PRSi) (7) 



e 



= —mt,{s,a^''Pnk)F^, (8) 



Stt^ 

= {^,ni,{s.<y"'WRh,)G%. (9) 



Here i and j are color indices (suppressed in color singlet terms), Prl = 
75)/2, and a^'^ = |[7'^,7'^]. The primed operators, which are not generated 
at leading order in the Standard Model, are obtained by taking L ^ R 
everywhere. Here we have ignored the electroweak penguin operators O7-10 
and the contributions to the dipole operators proportional to the s-quark 
mass, rris. 

Following the standard procedure for incorporating QCD corrections we 
match the Wilson coefficients at a high scale to loop diagrams containing 
heavy particles present in the full theory, and then use the renormalization 
group equations (RGE) to run the coefficients to the low scale where mesons 
decay. We incorporate leading order QCD corrections using the anomalous 
dimension matrices given in The initial conditions for the standard 
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Figure 1: Box and penguin contributions to the b sss transition. The 
bottom row shows contributions to the chromo-dipole operator. We show 
the mass insertions for pedagogical purposes but perform calculations in the 
mass eigenbasis. 

model coefficients are also in ^1]; only C2, and Cg are nonzero at leading 
order. Leading order running of the standard model coefficients has mixing 
between all eight operators, Ci_6,7,g, presence of a tree level 

contribution to 02- Since right-handed squark mixing only contributes to 
the primed operators and gives no tree level contributions, the leading order 
SUSY running is simpler: 0'^_q mix only amongst themselves, as do O'^ g.^ 

The SUSY contributions come from box, penguin, and dipole diagrams. 
Figure H shows sample diagrams with mass insertions schematically indicat- 
ing the mixing. However, since we are allowing for large mixing between the 
2nd and 3rd generation squarks, we use the mass eigenbasis for our compu- 
tations. Furthermore, we find that the region where the effect on S^k in the 

— i> (f)Ks channel is maximized is a region where the squarks are highly 
non-degenerate, again calling into question the validity of the mass insertion 
approximation. ^ 

^However, we have checked that even including mixing between these sets of operators, 
which is formally at higher orders in as, does not affect our numerical results significantly. 

^Indeed, our results differ somewhat from previous investigations done in the mass 
insertion approximation, e.g. |12j . 



4 



The squark mass matrix we consider is motivated by models where a large 
right-right mixing between the second and third generations is expected, such 
as in ,6J where this mixing is related to the large mixing in atmospheric 
neutrinos. In addition to this new contribution we must include other off- 
diagonal terms that already exist in the Minimal Supersymmetric Standard 
Model (MSSM), namely gnqi couplings induced by the cross term between 
the Yukawa couplings and the fi term.^ Of the down type squarks only the 
third generation can have appreciable left-right mixing, which is proportional 
to rribfi tan (3. Thus the mass matrix takes the approximate form 



/ mirf 

ml, 

m|rf 



\ m^/i* tan /5 m*^ 

We define the mass eigenvalues and mixing matrices as follows: 

U^MjU = diag(mii, mis, ^h, ^ri,^R2, ^m), (H) 

where U is a unitary rotation matrix. Without loss of generality we assume 
that m/j3 < m^2 ("we allow an arbitrary mixing angle). In the mass eigenbasis 
3x6 mixing matrices, F^'^, appear in the quark-squark-gluino vertices. They 
are related to U by 

UaA = ( r| ) , (12) 

where i = d,s,b labels the gauge eigenstates, A = 1, . . . , 6 labels the mass 
eigenstates, and the index a labels states in the basis {(11 sl bi sr 6r). To 
investigate the effect of 2nd and 3rd generation mixing of the right-handed 
squarks we parameterize the mixing matrix as follows: 

U = $(05, 06) RseiO^e) i?35(^35) ^^56(^56), (13) 

where $ = diag(l, 1, 1, 1, e*'^'^, e**^") is a phase matrix, and Rij{6ij) is a 2 x 2 
rotation in the ij plane. The angle 6*35 can be solved for using our assumption 

■^For simplicity we ignore terms that may arise from trilinear soft terms. 






mfoyutan/? 




m 



Rsb 



^2 



m 



Rb 



(10) 
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that there is no mixing between and sr. The Wilson coefficients in the 
mass eigenbasis were previously given in [TTl HH] . We reproduce them here, 
correcting two typographical errors. 



/ . ^ sA^ hA^ sB'- si 



B 



\AB 



1 5 

'-Bi{xa,Xb) - -B2{xa,xb) 

y y 



+ 2^ J- saJ- 



bA 



(14) 



a 



ml 



a 



/ . ^ sA^ hA^ sB'- sB 
\AB 



7 1 

'-Bi{xa, xb) + -B2{xa, Xb) 



^R*tR 



^sA^ bA 



1 3 
-Ci{xa) - 7;C2ixA) 

/ 



^R*-rR 



I sA^ bA^ sB'- sB 



y io 



^R*-rR 



-^C'i(a;A) + ^C2(xA) 



(15) 



(16) 



2 sA-"- bA'- sB'- sB 

m ' 



2 7 
--Si(xa, Xb) + 77S2(xa, Xb) 
o o 



1 3 

-Ci{xa) - -C2ixA) 
/ 



^R*-pR 

^^ b 

^R*-pR 



4 

--L)i(xa) 



^R*-pL 
{ J- 



^sA'- bA 



A 

m 



1 3 

■-Di{xa) + -D^{xa) 
2 



1 3 
'-D2{xa) + -/^4(a;A) 
o 2 



-g^2(a;A) 



(17) 
)(18) 



(19) 



Here we use the definition xa = ""^a/^S' where rrig is the gluino mass. The 
loop functions are given in Appendix El Note that the contributions due to 
left-right mixing only enter the dipole operators C'^ and C'^ where they are 
enhanced by a factor of rrig/mb over the right-right mixing contributions to 
the same operators. Also notice that with our choice of the mixing matrix, 
U, the mass eigenvalues rfiLi and rhm do not enter the Wilson coefficients. 
Because there was some disagreement in the literature, we have explicitly 
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recomputed the box contributions (proportional to Bi and i?2)- However, the 
penguin contributions (proportional to Ci^2 and -Di_4) are well established. 
See for example, [TJI. We found several inconsistencies in the literature which 
can be remedied as follows. In Equation (A. 8) of jll] the coefficient 
of Bi should be replaced by — and the expression for in Equation 
(A. 16) should be multiplied by — |. In Equation (42) of [12] the factor of ^ 
should be |, and in Equation (43) the factor of ^ should be |. The loop 
function M2{x) in ^H] should be multiplied by —x instead of — ^. Finally, 
in |H] there are typos in each line of Equations (B.4a-e) and in (B.6a,b). 

2.2 B^,^cl)Ks 

We now specialize our discussion to the B^ (pKs decay, with the goal of 
computing the contribution to the CP asymmetry measured in this chan- 
nel. In addition to the Wilson coefficients we have presented, we must also 
compute the hadronic matrix elements of the operators. The calculation of 
these matrix elements is non-perturbative, so approximations must be made 
in order for us to make progress. In the naive factorization approximation we 
break each matrix element up into a pair of color singlet currents, one which 
creates the (f) from the vacuum and the other that mediates the B^ ^ K 
decay, and we discard any color-octet currents. 

For the operators Os'lg there are two ways of contracting the external 
quarks with the quark fields in the operator. After employing Fierz trans- 
formations as necessary, and using the identity 5ij5ki = j^SuSkj + 2T°^T^j to 
rearrange color indices to form singlet currents, we arrive at the following 
matrix elements [S]: 

{<PKs\OsAB',) = \H(l + j^j (20) 

{<PKsmB^,) = \h (21) 

ms\0,m = ii/i- (22) 

where H = 2(e<^-pB)/(^m^F+(m|) (see Appendix iBlfor definitions of the decay 
constant and form factors). The same results hold for the matrix elements 
of the corresponding primed operators because the axial vector currents do 
not contribute, so the chirality of the operators is irrelevant. We take = 3 
throughout our analysis. 
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The matrix elements of the chromo-dipole operators O^'^ are more difficult 
to analyze, so we show the details explicitly in Appendix IHl following [S]. 
These manipulations yield 



(4'Ks\0,t\Bl) 



2-K 2m 



(23) 



where our definition of k, agrees with ^T] up to a sign convention.^ Nu- 



merically we find K = — I + (9 
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-1.1. However, there were many 



assumptions about the quark momenta that go into this estimation of k, so 
the numerical value of —1.1 should be taken as a guideline only. We will 
present our results for various values of k to demonstrate the dependence. 

There is one final ingredient in the Standard Model contribution to the 
amplitude. This comes from the one-loop matrix element of O2 when the 
charm quarks are closed into a loop. It is given by ^7] P 
G{mc, /i, g^)) with 



aa^-i (10 I 



dxx(l — x) In 



— x{l — x)q^ 

7 



(24) 



1.35 GeV and = ^ which gives P 



-0.015 



Numerically we use 
0.011i)C2. 

Putting the factorized matrix elements together the amplitude for Bd 
4>Ks can be written 



H 



+ H 



9 2tt 



2m 



^ ^ C 



(25) 



ISUSY 



H 



+ H 



1 1 



-A^ + J^](C', + C',) + -C', + --C', 



ANr 



1 



K 



2ti 2m 



(26) 



^Notc that in JT] the overall factor of j in the first line of Equation (18) should not 
be multiplying the last term that contains kdm- 
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The time-dependent CP-asymmetry is given by 



a^K{t) = C^K cos(AMBot) + S^K sin(AMBot), (27) 

where 



C -1^ S -1}^ (28) 



Here A is defined as 



The ratio q/p from mixing is dominated by the standard model and is 
nearly a pure phase, e*^^, where (3 is the standard angle of the unitarity 
triangle. In the Standard Model the ratio of amplitudes A/ A is real, i.e. 
there is no CP violation in the decay, rather all CP violation results from 
mixing. On the other hand, 0{1) phases in the supersymmetric contribution 
can give the ratio a phase, 0^^^^. Then we have 

X = ei(2/?+0su«-)M ^ s^K = sin(2/3 + ^^usy). (30) 
\A\ 

Thus the presence of a phase in the down squark mixing matrix can alter 
the measured value of S^k from the standard model prediction of S^f^K = 
sm2Pj/^K = 0.73. The amount of deviation is described in Section 1231 

Also note that P possesses a strong phase that is not present in 
The presence of a weak phase in ^^^^"^ then allows for the possibility of 
nonzero direct CP-violation, namely C^k 7^ 0. We do not pursue this sig- 
nature further here, as quantitative statements are difficult due to the large 
hadronic uncertainties. 



2.3 Constraints from 6 ^ 57 

A large mixing between right-handed strange and bottom squarks generates 
the operator O'^ through penguin diagrams, as in Figure |21 Therefore the 
tight experimental constraints on the branching ratio BR[B -^s7) serve 
to limit the contributions from squark mixing. 

In the model we consider there are two important contributions to the O'^ 
operator. We can classify the contributions according to where the helicity 
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Figure 2: Sample contributions to the 6 — 57 transition. We show the mass 
insertions here for pedagogical purposes, but perform calculations in the mass 
eigenbasis. 

flip for O'^ arises. In the first contribution, the helicity flip is present on the 
external 6-quark line, and gives a contribution proportional to the 6-mass. 
This contribution is present even when the only mixing between the squarks 
is an off-diagonal mixing between the right-handed squarks of the second 
and third generation. The constraint on this contribution is relatively mild.^ 
The second contribution has a helicity flip on the gluino line, so is enhanced 
relative to the first contribution by a factor of nig/miy (see Equations ()18p 
and ()19|)). This contribution is only present if there is left-right mixing in 
the squark matrix. Because of the rrig/mf, enhancement, this contribution 
is relatively strongly constrained. In our framework, this contribution arises 
only from the combination of a left-right mixing between the b squarks and 
the right-right mixing between the b and s squarks. The result is that for 
large values of /itan/5, a smaller bR-SR mixing is allowed. 

When there is no significant off-diagonal mixing among the left-handed 
squarks, we can write: BR{b 57) oc |C^P + IC'^l"^, where the first contri- 
bution is from the standard model and the second is from supersymmetric 
penguins.^ 



A recent theoretical evaluation within the Standard Model gives |18j 



After rescaling to limit the photon energies to > 1.6 GeV (for details see 
|18j). an averaging of experimental results from BaBar, BELLE, CLEO, and 

'^We note that the limits of reference assume that new contributions to 6 ^ 57 are 
summed incoherently. In general, this will underestimate the contribution of the LL mass 
insertion. 

^The two contributions are added incoherently because they contribute to different final 
helicity states of the s-quark. 




BR{b s-f)th = (3.60 ± .30) x 10 



,-4 



(31) 
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ALEPH [in] yields 



BR{b s-f)exp = (3.29 ± .34) x 10"^ (32) 

The experimentally measured branching ratio is actually shghtly smaller than 
the standard model prediction, which leaves little room for new physics con- 
tributions. Subtracting experiment from theory we find: 

BR{b s-f)th - BR{h ^ S7)e.p = (.31 ± .45) x 10"^ (33) 

We will require that the supersymmetric contribution keep the theoretical 
prediction within 2a of the experimentally measured value. This means ad- 
ditional contributions from supersymmetry can be roughly 1/6 of those in 
the Standard Model. For simplicity, and to avoid the theoretical uncertainty 
associated with the direct calculation of the branching ratio, we will con- 
strain the supersymmetric contributions by requiring \C'^\'^ < 0.16 x |C^p 
where both coefficients are calculated to leading order. Thus we are making 
the simplifying assumption that the higher order QCD corrections affect the 
two operators in the same way. 

2.4 Numerical Analysis 

Within the framework we have chosen, motivated by atmospheric neutrino 
oscillations, there are four mass eigenvalues, two mixing angles, and two 
phases in the down-squark mass matrix that enter in the computations of 
S^K and BR{b 57). However, the fact that the neutron electric dipole 
moment (EDM) has not been observed strongly constrains the phase of 
especially for large fi tan (3. We have checked that allowing a non-zero phase 
of fi does not substantially affect even our quantitative conclusions. Therefore 
we conservatively take the phase of fi to be zero for the remainder of this 
paper. Including the gluino mass we are then left with eight essentially 
unknown parameters. In order to reduce the size of the parameter space 
we will investigate two limiting cases: the case where the (flavor-diagonal) 
mixing between left- and right-handed b squarks is negligible and the case 
where such mixing, when coupled with the large right-right mixing, leads to 
the dominant contribution. We refer to the latter CclSG clS ^^LR + RR mixing." 

Note that by taking /i to be real, the remaining phase only appears as an 
overall phase in Tji. As a result, all SUSY diagrams have the same phase. 
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2.4.1 Dominant Right-Right Mixing 

First we consider the situation where the contribution from bi-bn mixing is 
negligible, i.e. mf,/i tan /5 <^ m^, the mass scale of the squarks.'' Here the 
parameter space is reduced: in this limit there is only one mixing angle, 
one phase, and three mass eigenvalues that enter the computation of the 
Wilson coefficients. The presence of mixing with an order one phase in the 
right-handed down squark sector can significantly alter the measured value 
of S^K- Our first question is whether a large right-right mixing between the 
down squarks can reproduce the central value for sin 2f3 in the (pKs channel 
measured at the 5-factories. We find that using the central value of our 
estimate for k = —1.1, it is possible to reproduce the observed central value 
and accommodate the constraints from b S7. However, this estimate 
for K, is highly uncertain, and increasing the magnitude of k increases the 
contribution to S,j,K without changing the contribution to 6 ^ 57. Therefore, 
we present our numerical results for two 1.1, and a value with 

greater magnitude, k = —2, which we still view as reasonable given the 
substantial uncertainties involved in its estimation. 

In Figure 01 we show contours of S^k as a function of the gluino mass 
and m/j3. We have also chosen values for the mixing angle and phase in 
which give the greatest deviation of S^k from the Standard Model prediction. 
Also shown are contours of the percent increase in BR{b 57) due to new 
physics and the corresponding values of AMbs (the latter will be discussed 
in Section El). For gluino masses around 200 GeV, S^k can take on values 
as low as —0.4 for n = —1.1 while still keeping the increase in BR{b S7) 
below 16%. For k = —2.0 the value of S^px can reach all the way to —1.0. 

Generally speaking, lighter squark and gluino masses increase the effect 
of the new physics contributions, allowing S^i^k to depart from the standard 
model expectation. But at the same time this increases the contribution to 
BR{b — i> 57) and runs up against the experimental constraint. 

In Figure |3] we plot the same contours as a function of the gluino mass 
and the heavier squark mass m_R2; with m^js = 300 GeV. For k = —1.1 there 
is a range of gluino masses where S^k can be below zero while still satisfying 
the 6 — > S7 constraint. The minimum possible S(I)K allowed by this constraint 

^Within the MSSM /xtan/3 cannot go to zero while satisfying experimental constraints. 
We find, however, that there is a portion of parameter space above the smallest exper- 
imentally allowed value of /itan/3 where the right-right mixing diagrams are dominant. 
Furthermore, right-right mixing dominates when rfiR <C ttt-l. 
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Figure 3: Contours of S^k (thin blue), percent increase in BR{h S7) (thick 
green), and AMb^ (dashed red) for Thfi2 = 5 TeV, optimal mixing angle and 
phase, and k = —1.1 (left) and k = —2.0 (right). 

decreases as the magnitude of k increases. For example, with k = —2.0 the 
minimum value of S^f^K is roughly —0.8. 

Finally in Figure El we plot the same contours as a function of the gluino 
mass and k, with fhff2 = 5 TeV and rfiRs = 300 GeV. The b ^ s'y constraint 
is independent of k, and we see how S^px decreases with the increasing mag- 
nitude of K. 

2.4.2 Dominant LR + RR Mixing 

Now we consider a second limiting case, where the diagonal left-right mixing 
leads to the dominant contribution to the 6-quark decays. The contribution 
from LR+RR mixing is enhanced in the dipole operators, so we may focus our 
attention on the coefficients Cg and C'^. To evade the constraint from b 57, 
while simultaneously getting a large effect in S'^i^-, we want to minimize the 
ratio C'^/C'g, which can be done by taking larger values of x, i.e. squark 
masses much heavier than the gluino. 

In Figure ini we reproduce Figure El but with a large value of /itan/3 ~ 35 
TeV. In this region the lightest squark mass eigenvalue, rhus, needs to be 
above 1 TeV to avoid the bound from b — > 57. For k = —1.1 the smallest 
S^K can be is about —0.4 for a small gluino mass, while for k = —2.0 any 
value is possible. In this case the ^-dependence is very simple because the 
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Figure 4: Contours of S^k (thin blue), percent increase in BR{b 57) (thick 
green), and AMb^ (dashed red) for rhus = 300 GeV, optimal mixing angle 
and phase, and k = —1.1 (left) and k = —2.0 (right). 
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K 

Figure 5: Contours of S^k (thin blue), percent increase in BR{b — > 57) 
(thick green), and AMb^ (dashed red) for fhR2 = 5 TeV, = 300 GeV, 
and optimal mixing angle and phase. 
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Figure 6: Contours of Sfj,K (thin blue), percent increase in BR{h 57) (thick 
green), and AM^^ (dashed red) for 177.^25 ""^ls? ^R2 ~ 2 TeV, optimal mixing 
angles and phases, /itan/3 ~ 35 TeV, and k = —1.1 (left) and k = —2.0 
(right). 

main contribution to S^k comes from a single operator O'g whose contribution 
is directly proportional to k. Thus an increased absolute value of k directly 
increases the effect in S^j^K without affecting the bound from h 57. 

2.4.3 Combination of Contributions 

In order to ascertain the relevance of each of these two regimes we scanned 
over the parameter space searching for the minimal values of S^px as a func- 
tion of the product /itan/5. The result is shown in Figure[Z|for sevaral values 
of K. From the figure it is clear that for k = —1.1 there is a slightly larger 
effect on S^k for larger values of /itan/5, though the entire region allows for 
a substantial deviation from the standard model. Notice that once k reaches 
—2.0 any value of S^k is possible. 

To get a sense for the relative size of the RR and LR + RR contributions 
to S^K, we can compare the magnitude of the two terms comprising Cg, 
Equation (jl9|) . Not surprisingly, for /itan/3 greater than about 25 TeV the 
LR + RR contribution dominates by an order of magnitude. However, even 
for /i tan /3 < 1 TeV there can be points where the LR + RR contribution 
to the chromo-dipole operator is just as important as the RR contribution. 
This underscores the importance of treating this calculation in the mass 
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Figure 7: Minimum values of Sij,K as a function of yU tan /3 for various values 
of n resulting from a scan of 58 million points that satisfied the h ^ 
constraint. Note that n = —2.0 can give Si^k = —1.0 for all values of /xtan/?. 

eigenbasis. 



3 Bs Mixing 

Mixing between b^i and sr also leads to a significant contribution to Bg-Bg 
mixing. In our scenario the effective Hamiltonian that receives such contri- 
butions consists of three operators that have nonzero coefficients: 

(34) 

The Wilson coefficients at the high scale are obtained by matching the ef- 
fective Hamiltonian to the AB = 2 squark-gluino box diagrams like those 
shown in Figure |H1 The result is given by^: 

2 



2 / y s ^ 'I'l 
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LL 



"^9 AB 



lA-"- bA^ sB^ bB ^ 



-B2{xa,xb) 

O 



(35) 
(36) 
(37) 



where the loop-functions are defined in Appendix ^ the F's are defined in 
Equation (HU, and again xa = m\/rn?~. 

^This result differs from 20 , but agrees with subsequent analyses, e.g. 
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Figure 8: Diagrams contributing to i?s-niixing through large bR-SR mixing. 
The diagram on the right induces effective bi-SR mixing which contributes 
toCf^andCf^ 



The leading Standard Model contribution to Bg mixing is induced by a 
top quark box diagram which yields the following effective Hamiltonian 



n'J' = cil^isrPMsi.PLb), (38) 

with the Wilson coefficient matched at Mw 

CIl = {V*V,,f S{ml/Ml), (39) 



167r2 



where 



^("> = 4(1 -xp — 2(r^- 

Before taking the hadronic matrix element of the effective Hamiltonian we 
must first take QCD corrections into account by using the renormalization 
group equations (RGE) to evolve the Wilson coefficients down to the low 
scale. The general NLO running of the Wilson coefficients for a AB = 2 
effective Hamiltonian is given in [21] and [221 and involves mixing among 
different coefficients. In our case only the two scalar left-left operators mix, 
while the vector right-right and the vector-left-left coefficients simply scale 
multiplicatively. For simplicity we have evolved both operators from Mw 
down to the mass of the 6-quark.^ 



^The supersymmetric contribution should in fact run from ^ rrig down to mf,. In this 
approximation we are ignoring corrections of order 1 — as{mg)/as{Mw) and potential 
contributions from the top quark in loops which is smaller. These corrections are part of 
the systematic uncertainty in our calculation. 
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The hadronic matrix element of the effective Hamiltonian between Bg 
and Bg states was calculated on the lattice [2S1 

{Bg\{s^'^Pnb)is^,PRb)\Bg) = ^msJlBM (41) 

{Bg\isPLb){sPLb)\Bg) = -7^(^^) nisJlBM (42) 

24 Vmf, + iTLsJ 

1 / M \ ^ 

m{s,PLb,){-s,PLh)\Bg) = ^[;;;;^J ^bJIBM, (43) 

with Jb, as given in Tabled Bi{mb) = 0.87(4)1^ ^2(7715) = 0.82(3)(4), 
and B^lmh) = 1.02(6)(9), where the first error is statistical and the second 
is systematic, excluding uncertainty due to quenching. The quark masses 
in the above expression should be evaluated at the scale /i. The hadronic 
matrix element for the left-left current operator of the Standard Model is 
identical to that of the right-right operator shown in Equation (j4H) . Finally, 
we can write the expression for the mass difference between Bg and Bg as 



AM 



Bs 



{Bg Hf/=2 Bg) . (44) 



AB=2 



The Standard Model and supersymmetric contributions interfere, 

The input parameters used in the calculation are given in Table d Our 
results should be compared to the standard model prediction which can be 
obtained roughly by taking Ti.^^"^ = Tifff^ in Equation (jUj) . which yields 
AMJ^^ ~ 13.8 ps^^. A more rigorous treatment given in |21j yields 

AM|f = 17.3lli Ps"'- (45) 

However, given the substantial uncertainty in the lattice evaluation of, e.g., 
Jb^ , it is probably appropriate to inflate this error, likely to the 25% level |2H] . 
The current experimental limit, combining results from the LEP experiments 
and SLD, is [25] 

AMb, > UA ps"^ (95 % confidence level). (46) 

Current and upcoming experiments are expected to be sensitive to mass dif- 
ferences much greater than the Standard Model prediction shown in Equa- 
tion ()45|) . At Run II of the Tevatron fI7] CDF is expected to probe up to 
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Parameter 


Value 


Parameter 


Value 


Parameter 


Value 


rrib 


4.2 GeV 




5.379 GeV 




-0.04 




174 GeV 




204 MeV 


as{Mz) 


0.1185 


Mw 


80.4 GeV 


Mz 


91.2 GeV 


^Bs 


1.461 ps 



Table 1: Input parameters used in the calculation of and AM^^. 



AM^^ of 41 ps^^ while BTeV is expected to achieve sensitivity to values 
up to AMs^ ~ 55 ps^^ Any evidence that AM^^ > 25 ps~^ from these 
experiments would be a clear signal of new physics. 

To illustrate our results we add contours of constant AM^^ (red dashed 
lines) to Figures EHHl In the case of dominant right-right contributions, i.e. 
small /i tan /5, the trend is similar to that of the previous section; lighter gluino 
and squarks give a larger SUSY contribution and thus increase AM^^ . Note 
that the supersymmetric contribution to the mass difference dominates over 
the Standard Model in significant regions of the supersymmetric parameter 
space, easily allowing AM^^ > 100 ps~^ where right-right mixing dominates. 
Such values for AM^^ are certainly beyond the reach of the experiments 
mentioned above. 

In the case of dominant LR + RR mixing the modification of mixing 
is not as striking. In the example given in Figure IHl values of AM^^ are much 
closer to the standard model prediction. Restricting ourselves to areas that 
respect the 6 — *• S7 bound gives a yet lower value, within the reach of upcom- 
ing experiments. We should point out that this is not generic since Figure El 
only represents a slice of parameter space. Other choices of parameters can 
give higher values of AM^^ (above 30 ps~^) for high values of /itan/3. The 
correlation of these results to those in (t>Ks will be discussed in the 

next section. 

Finally, we should comment about possible CP violation in the Bg-Bg 
system. In the standard CKM scenario the Bs Bg amplitude does not have 
a CP violating phase (in the Wolfenstein parameterization), so no indirect 
CP violation is expected. In our scenario, however, Bg mixing can involve 
the phases from the down-squark mass matrix. In the cases where the SUSY 
contribution to Bg mixing dominates the SM, measurements of CP violation 
in Bg —>■ J/ipcf) will be sensitive to these phases. 
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4 Correlation 



In this section we will discuss the correlation between S'^x and AM^^ in 
the context of large b-s mixing. Because the effect on Bs mixing is very 
different for the two limiting regions of parameter space, we will discuss 
them separately. For related studies see [T21I2H1- 

4.1 Dominant RR Mixing 

In this region of parameter space the operators 0'^_q make large contributions 
to Sfpx, while there is essentially only one contribution to Bg mixing, namely 
that from the operator shown in Equation (j35|) . Unfortunately there is no 
simple, precise, relationship between the combination of the AB = 1 oper- 
ators and the operator responsible for Bs mixing. In general, they depend 
quite differently on loop functions. 

In spite of this, one can make the following strong statement. In cases 
where there is a large shift in S^k away from the Standard Model expectation 
due to the operators 0'^_q, and the RR contribution to O'g (the dominant 
right-right mixing scenario), there is a large contribution to Bg mixing. To 
see this, we first note that the squarks and gluino must not be too heavy, and 
the b-s mixing must be large in order to have a large contribution to S^k- 
This suggests a minimum contribution to the Bg mixing. However, there 
is the worry that it might be possible to fine-tune parameters to somehow 
drastically suppress the contribution to Bg mixing; for example, by choosing 
squark and gluino mass ratios to minimize the value of the functions Bi and 
B2 in Equation (j35p . We find that this is not possible, however. In order 
to have a very large contribution to S'^x, one is pushed into a region of pa- 
rameter space where the gluino, and at least one of the down-type squarks is 
light. Furthermore, the splitting between this light squark, which represents 
a mixture of b and s squarks, and the masses of the heavier squarks must be 
large to avoid a super-GIM cancellation. Once this qualitative picture for the 
spectrum is identified, it is easy to check that there cannot be a cancellation 
of the contribution to Bg mixing in this case. 

Now we present our results quantitatively. Since our goal here will be 
to show that large deviations in S^k will correspond to large contributions 
to Bg mixing, we plot the minimum achievable value of AMb^ for a given 
value of SfpK- The minimum is found by scanning a parameter space that 
consists of the parameters {m^, 771^2, ''^H2, '"^ra, cos 6*56, 06 }• As discussed in 
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Equation (jl2j) . cos represents the mixing angle between the right-handed b 
and s squarks, and 06 represents the phase corresponding to this off-diagonal 
term. As a parameter space, we take:^° 





e 


(200, 700) GeV 


rhL2 


e 


(300, 2500) GeV 


rhR2 


G 


(300,2500) GeV 




G 


(250, 1000) GeV 


(^56 


G 


(0,77/2) 


06 


G 


(0,27r). 



The lower limits on the masses are motivated by direct searches, while the 
upper limits are motivated by naturalness considerations. A scan would 
generate a scatter plot of S^k vs. AMb^- For a given resultant value of S'^x, 
we find the combination of parameters that yields the smallest contribution 
to Bs mixing. This is essentially equivalent to taking the boundary of the 
region generated by the scatter plot. 

As discussed in Section |21 there is considerable dependence on the variable 
K, which has a relatively large uncertainty. So we repeat the above exercise for 
several values of k, displaying the results in Figure IHl Adding the constraint 
from b ^ s'y modifies these results as shown in Figure fTUl The contours in 
FigurelTUlnotablv do not extend as low in S^k because the b ^ sj constraint 
removes the region of parameter space that allowed us to obtain those values 
in Figure ini 

The take-home message from the figures is a simple one. If the hint of the 
deviation in S^k measured in the B — >• (f)Ks persists and it is attributable 
to a scenario with dominant RR squark mixing, it will result in a large 
contribution to Bg mixing, which will be a clear indication of new physics 
observable at the Tevatron. 

4.2 Dominant LR + RR Mixing 

In the region of parameter space where /itan/5 is relatively large, the expec- 
tation for AMb^ is very different. In this region the main contribution to 
S'^x comes from the LR + RR contribution to the dipole operator Og. This 

"'^^Here we have taken /i tan = for simplicity. Nonzero values can weaken the corre- 
lation somewhat, as shown in Figure ITTI 
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Figure 9: The minimum value of AM^^ for a given value of S^k- It is 
found by scanning over the parameter space given in Equation (j47p . The 
corresponding curve is shown for several values of k, the coefficient of the 
dipole operator, as defined in Equation ((221) • 
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Figure 10: The same as Figure El but with the 6 — 57 constraint applied as 
discussed in the text. The minimum value of AM^^ for a given value of S^k- 
It is found by scanning over the parameter space given in Equation ()47|). 
The corresponding curve is shown for several values of the coefficient of 
the dipole operator, as defined in Equation (|23|l . 
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contribution can be sizeable even when the squarks and gluinos are heavy 
(squarks can be at the TeV level or higher). This is the significant difference 
between the two limiting cases. Heavy squarks and gluino mean that the 
contributions from the operators 0'^_q are small. Similarly, the operators re- 
sponsible for Bs mixing, which come from box diagrams and resemble 0'^_q, 
can also be small. The bottom line is that a large contribution to S^j^x is 
possible without a large addition to AM^^. This is borne out numerically, 
as shown in Figure where the points allowed by the 6 — > 57 constraint all 
give AMb^ very close to the Standard Model expectation. 

We have seen that the LR + RR contribution to SfpK can be important 
even for fairly small values of /itan/5, so it is natural to wonder what con- 
clusions can be drawn about AM^^ in the regions where both contributions 
are important. To answer this question we again performed a scan of the pa- 
rameter space, this time collecting points with maximal and minimal AMb^ 
as a function of /itan/5 with the additional requirement that S^k < —0.2 
for the nominal value k, = —1.1. The results are shown in Figure ITTl In ac- 
cord with what was stated above, points with the largest /itan/3 give smaller 
contributions to Bg mixing. In fact, for /itan/3 > 40 TeV any effects on 
AMbs will be indistinguishable from the Standard Model expectation. This 
apparent upper bound may be interpreted as follows. For large /xtan/3 the 
severe 6 — 57 bound is pushing us to regions where RR mixing is small or 
masses are high, both of which disfavor large contributions to Bg mixing. At 
the other end of the spectrum, for /itan/5 < 2.4 TeV all points in the scan 
gave values of AMb^ > 30 ps~^, a clear signal of new physics above even 
the largest Standard Model predictions. This trend continuously connects us 
back to the result of the previous subsection. 

Mixing between the first two generations of squarks has fallen outside the 
main scope of this paper. We should mention, however, that larger values 
of the lightest squark masses, around 500 GeV, may well be preferred by 
constraints from K-K mixing. To see this, note that if we believe that the 
Cabibbo angle originates through the down Yukawa matrices, then the down 
Yukawa matrix has the structure: 



where h is some unknown Yukawa coupling and Ac is the sine of the Cabibbo 
angle. Diagonalizing this matrix requires a rotation on dL and sl of 0{Xc) 




(48) 
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Figure 11: Maximum and minimum values of AMb^ as a function of ^tan/? 
resulting from a scan of 150000 points in the parameter space with S^k < 
—0.2 and k = —1.1. The plot on the left shows an enlargement of the region 
of small /X tan /3. 



which induces an off-diagonal element of order h^Xc- Using the phenomeno- 
logical relationship, md/ms ~ Xq, we find that a rotation between (Ir and sr 
of O(A^) is needed to complete the diagonalization. Then, due to the lack of 
degeneracy between the squarks, the induced d-s mixing can lead to a large 
contribution to the K-K mixing. This suggests that heavier squark masses, 
perhaps above 500 GeV are preferred, barring some accidental cancellation 
with the (unknown) (1, 2) element, h, in the above matrix. 

By the above reasoning, if one wishes to achieve an S'^x that differs 
significantly from the value of sin 2f3 as measured in the in the J / ipKs 

channel, there may a theoretical prejudice to prefer scenarios where large 
squark masses are more easily accommodated, such as the LR+RR dominant 
case or the RR dominant case in conjunction with a large value of k. 



5 Conclusion 

There exist a class of models, motivated by Grand Unified theories and the 
large observed mixing in atmospheric neutrinos, where it is natural to have 
a large mixing between the right-handed b squark and s squark. We have 
found that there exists a range of parameters where such mixing induces a 
significant deviation in S^k from the Standard Model expectation of sin 2/3 
as measured in the channel 5° J/ipKs- 

In particular, the central value for Sif^K from BaBar and BELLE can be 
accommodated without conflicting with the measured value of fe — > 57, using 
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the naive estimate for k, the hadronic matrix element for the chromo-dipole 
operator. For larger values of k, any value of S^px is allowed within the 
constraint. 

There are two possible origins of a substantial modification to S^k- The 
first solely involves a large right-right mixing, with no contribution from the 
mixing proportional to rribfi tan /5. In this small gluino mass is required, 

near the experimental bound. Correspondingly, there is a large contribution 
to Bs mixing, a consequence which will be testable at the Tevatron Run II. 
In the second case, we consider the mixing from the combination of the large 
right-right mixing and a large rubH tan (3. this case, squarks and gluinos 
need not be light, so Bg mixing need not be large. In particular, for very large 
values of /itan/3 the prediction for Bg mixing is indistinguishable from the 
Standard Model prediction, when current errors on lattice matrix elements 
are taken into account. However, a substantial improvement in the Standard 
Model prediction for Bg mixing still may allow an effect to be seen at the 
Tevatron in this case. 

Note Added 

While completing this paper we received References fl^ I30j. There is some 
overlap with these papers, which also consider supersymmetric contributions 
to 5° ^ (l)Ks. 

Regarding j^H], in places that we overlap, we agree qualitatively with their 
results, though there may be some quantitative differences. These are likely 
due to the fact that they work in the mass insertion approximation. Indeed, 
allowing large mixings and hierarchies that cannot be described by mass 
insertions gives us larger contributions in the 'pure' RR mixing case. Other 
possible differences may arise from a different treatment of the hadronic 
matrix elements and the fact that constraints from b s'j were not imposed 
in the same way. 

We differ from both papers in our emphasis on the mixing induced by a 
combination of mhfi tan j3 along with a large flavor-changing RR element in 
the squark mass matrix. Though double mass insertions are briefly discussed 
in [21], both [221 and [30] focus on the contributions of single mass insertions, 
including flavor off-diagonal LR mixing. 

We would like to emphasize that the double mass insertion {LR + RR 
mixing) does not necessarily describe the same physics as a single LR flavor 
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mixing insertion, a point also mentioned in • Treating a LR + RR mixing 
as a pure may miss important contributions due to the RR mixing only. 
For example, we find maximal values of AMb, for intermediate values of 
H tan (3 which are much higher than those plotted in j^Hl with pure LR mixing. 
This might be due to a RR mixing contribution which is sub-dominant in S^k 
but which nevertheless gives a large contribution to AM^^. This example 
illustrates how our framework differs from analyses which consider only one 
mass insertion at a time. 
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A Loop Functions 

We include the loop functions for completeness. We use the same definitions 
as in Here xa = ''^|^/'"^|- 




x\ log Xa x% log xb 



4:{xa- xb){xa- ly 4:{xb - xa){xb - ly 
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B2{xa,xb) = 
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Xa log Xa Xb log xb 
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B Chromo-dipole Matrix Element for Bjj 

(t^Ks 

In this appendix we show explicitly the computation of the matrix element of 
OW in the naive factorization approximation. The computation of a similar 
quantity for tt+tt" decay can be found in pT| . 

We start with 

= j^,m,{-s,a^^^Pnh)G% (57) 

and then connect a quark current through a virtual gluon to form a four- 
quark operator. This step depends on the convention used for the covariant 
derivative. Our convention is that = d^ + igT°'A'^^, and we have checked 
that this is consistent with the Wilson coefficients for both the standard 
model and SUSY contributions. This yields the operator 

^.4 = i-^^m,{-s,a^^^Pnh,){qul,T^fli) (58) 

where k = pb — Ps is the gluon momentum. In the naive factorization ap- 
proximation the color-octet current, {qk'J^iTl^eqe), cannot produce a physical 
0, so the (p must be produced by the s and q operators. To factor the matrix 

-'^-'^In particular, note that (351 appears to use the opposite convention of |31[I33| . 
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element we first use the equations of motion to simplify the tensor current. 
This yields: 



{qkl,T;:,qd (59) 



where the last term was simplified using the conservation of the quark cur- 
rent, k'^iqk'Jfj.Tj^eQe) = as in Then by a Fierz transformation and judi- 
cious use of Dirac matrix identities, this can be brought to the form [HIES EH] 

X [2mi,{sa^PLqi){qkl,,PLhj) - 'imhisiPRqe){qkPLbj) 
+2ms{sa''PRqi){qk'J^,PRbj) - Ams{siPLqi){qkPRhj) 
+2{p,)^{{siYPLqd{qkPRhj) + {siPRq^Wkl^PRh) 

+i{s,a>'''Pnqe){qkYPRbj) - i{sauPLqi){qkCJ^'' PRhj)}] 

Next we use the following parameterization for matrix elements. 

(0(p^,e^)|^7^s|O) = Umy^ (61) 
(0(p<^,e^)|ss|O) = (62) 

2777 

(0(p^,e<^)|.-a'^^.|O) = -^/^mj^(e^p;-e;p^). (63) 

Also, 

{K\pK)\-sYh\B'i{pB)) = {vb+VkYF+H) 

+ {PB-PKYF_{t) (64) 
{K\pK)\-sh\Bl{pB)) = 6^,'ABKF4t) 

+ 6^,\pB-pKrF4t) (65) 
{K'{pK)\sa'^'{l±l,)bmpB)) = -2zs{p^bPk - PbP'k) 

± Se'"'^''{pBXPKa-pBaPKx). (66) 

Here t = {pb — PkY, ^bk = it^'b ~ ^"k^ ^bs = mb — rus- Notice we have 
corrected the sign in Equation ()63|) compared to the similar expression in . 
Heavy quark effective theory gives the relation s = (F+ — F_)/Amh [SI]. We 



28 



also make the kinematic assumptions that the 6-quark carries all of the B- 
meson momentum and that the momentum is equally divided between its 
two constituent s-quarks. Thus pb = Pb and fc^ = — + 

Putting all the pieces together gives: 



+ 



ml 
2A;2 



"2P 



1 + 



m 



B 



mj( + m 



2 F- 



m 



K 



— 2m 




(67) 



Because the matrix element is nonsingular we have F^{0) = 0. Then for 



small t, due to simple pole dominance we have F_ 



F_(0) = 0. jSS] 



By ignoring the difference between 6-quark and 5-meson masses we arrive at 

the estimate cited earlier, k = — | + (9 ( —r I — —1.1. The sign convention 

differs from jH], and the slight difference in magnitude can be traced to our 
replacement of (pb + Ps)^ with 2p^ using the conservation of the added quark 
current. In IJJI a similar quantity S^j^K = |k is used. However, they quote a 
value S^K — —0.76 which appears to match iS^^r found in [^1]. This would 
correspond to a value of k ^ —1.71. 
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